EEE 


D 8 
Sy 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
Ll 
LLL IIIIIIIII BBBBBBBBBBBB RRRRRRRRRRRR TTTTTTTTTTTTTTT «LLL Ll 
LLL HII BBBBBBBBBBBB RRRRRRRRRRRR TTTTTTTTTTTTITTT «LLL + 
LLL LITT BBB RRRRRRRRRRRR TTTTTTTTTITITITT «LLL L 
LLL 111 BBB BBB RRR RRR TTT LLL | Ll 
LLL II] BB BBB RRR RRR TTT LLL LI 
LLL 11] BBB BBB RRR RRR TTT LLL | Ll 
LLL II] BBB BBB RRR RRR TTT LLL | Ll 
LLL Il] BBB RRR RRR TTT LLL Ll 
LLL Il] BBB BBB RRR RRR TTT LLL Ll 
LLL III B RRRRRRRRRRRR TTT LLL Ll 
LLL Ill BBBBBBBBBBBB RRRRRRRRRRRR TTT LLL LI 
LLL Ill 8 RRRRRRRRRRRR TTT LLL ul 
LLL III BRA BBB RRR RRR TTT LLL 
LLL II] BBB RRR RRR TTT LLL Ll 
LLL II] BER BBB RRR RRR TTT LLL Ll 
LLL II] BRB BBB RRR RRR TTT LLL Ll 
LLL III ERB BBB RRR RRR TTT LLL Ll 
LLL Ill BSB BB RRR RRR TTT LLL Ll 
LLELLLLLLLLLLLLL IIIIIIIII BBBBBBBBBBBB RRR RRR TTT LELLLLLLLLLLLLL Ll 
LLELLLLLLLLLLLLL TITIII111 BBBBBBBBBBBB RRR RRR TTT LI 
LLELLLLLLLLLLLLL LIIIIIII BBBBBBBBBBBB RRR RRR TTT LELLLLLLLLLLLLL 


| 
LLLLLLLLLLLLLLL | 
| 
| 
| 


**FILE®*1D**OTSPKDIVS 


000000 TTTTTITITT SSSSSSSS PPPPPPPP KK KK DDDDDDDD IIIIII VV 
000000 TTTTTTTITT SSSSSSSS_ PPPPPPPP KK KK DDDDDDDD III] VV 

00 00 TT SS PP PP KK KK ODD DD I] VV 

00 00 TT SS PP PP KK KK DD DD II VV 

00 00 TT SS PP PP KK KK DD DD I] VV 

00 00 TT SS PP PP KK KK DD DD I] VV 

00 00 TT SSSSSS PPPPPPPP KKKKKK DD DD I] VV 

00 00 TT SSSSSS PPPPPPPP KKKKKK DD DD II WV 

00 00 TT SS PP KK KK DD DD I] VV 

00 00 TT os PP KK KK DD DD I] VV 

00 00 TT 33 «6PP KK KK DD DD I] VV VV 

00 00 TT SS PP KK KK DD DD I] VW VV 

TT SSSSSSSS PP KK KK DDDDDDDD IIIII1 VV 

000000 TT SSSSSSSS PP KK KK DDDDDDDD III1I1 VV 

LL IIIT] SSSSSSSS 

LL IIIT] SSSSSSSS 

LL II SS 

LL I] SS 

LL I] SS 

LL I] SS 

LL I] SSSSSS 

LL I] SSSSSS 

LL I] SS 

LL I] SS 

LL I] SS 

LL I] SS 

LLELLLLLLLL III] SSSSSSSS 

LLLLLLLLLL HII! SSSSSSSS 


04-00 Page 
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OTSSDIV_PKSHORT 16-SEP-1984 00:31:31 YAX/VMS Macro V 
11801 ¢- 3 99:73:33 LIBRTL.SRCJOTSPKDIVS.MAR; 1 


SEP=1984 


-title osha iy phohert 
eident /1-001/ ; Edit 0G1001 


DADA DAARRRRAARRAAASAZASLLESELELLSE ALAS ELAR ARSE ESE REESE SER EERE SER SES ESSEC EES SY 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 


* 
* 
® 
® 
« 
* 
. 
* 
‘ 
* 
TRANSFERRED. a 
* 
* 
‘ 
* 
’ 
* 
« 
* 
* 
* 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
eaRPORAT Lon NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO ep 


IBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS 


NOT SUPPLIED BY DIGITAL. 


De Se Be Oe Se Se Se Se Be Oe Be Se Se Se Se Se Se Se Se Bese 
@eeeeneeeeneneee eee neenanne 


LARRARASASAASASASLAALASALALSSASAL ASE LES ARES SRLS RRS E SESE RRR ER RARER EER RRA SES SD 


a 
+ 


routine: 
OTS$DIV_PKSHORT 


facility: 
VAX/VMS OTS runtime Library. 
abstract: 


Runtime routine performs fixed decimal (packed decimal) division. 
The routine is called when precision and scale requirements for 
the quotient imply eulttere precision division. The routine is 
only called when such multiple precision division is required and 
when the divisor has a precision of less than 50 decimal digits. 
(Call ots$div_pk_long if mu!tiple preg ioten division is. 

required and the divisor has precision 30 or 31 decimal digits). 


author: Peter Baum 20-jun-1980 
modifications: 


1-001 Debess Grabazs 5 March 1984 
Made PLI routine into OTS routine. 


NAW WN OC OD NA UE WIN O OD NAME WIN O OO NAU EWR OO OO NOAU EW OOONOUS Wr 


DUPVDPWSVSVSISIIS BB BS BB EE ANNI ANIA PPI NINPONPIPONIPIAY 2 OO Os ee 


Be Se Se Se Se Ge Be Se Se Se Ge Se Se Se Se Ge Ge Se Ge Ge Ge Se Ge Ge Ge Ge Se Ge Se Ge Se Se 


OTS$DIV_PKSHORT 
1-601” 


P=-1984 00:31:31 VAX/VMS Macro V04-00 
eats 90:31:33 YMTGRTL. SRCIOTSPKDIVS MAR: 1 


documentation file: THEORY.MEM 
functional description: 


WR OOO NAUES WN —-OOONAUES WOO 


This routine calculates: 


z=x/y 

let a = scale(z) + scale(y) = scale(x) - 31 + grastas 
b = scale(z) + scale(y) = scale(x) + prec(x 
¢ = 31 = prec(x) 
d = 31 = prec(y) 


this routine is called if b > 31 and d > 1 


Prior to the call: 
f c not 0 then shift x left by c. 
Thus x is a 31 digit packed decimal. 


input: 
O(ap) # of arguments 
4(ap) address of dividend (shifted left by c) 
8(ap) address of divisor 
12(ap) precision of divisor (high order bytes zeroed) 
16(ap) address of quotient 
20(ap) precision of quotient (high order bytes zeroed) 
24(ap) a as defined above(high order bytes zeroed) 
28(ap) das defined above(high order bytes zeroed) 
output: 


quotient returned at address specified by 16(ap) 


variable usage: 


Pe ee tr eet Re a tet Ot tet try tet eet ER PRN RE OS Or RN RA 


Se Se Se Se Se Se Ge Ge Ge Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Be Ge Ge Ge Ge Se Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Be Se Se Se Se Se Se Se Se Se Ge Sse Sse Ge Sete Se 


PWD OC OODNOAU EWN O OONOAUE WN - OOONOAM 


a bs tt a 


size 
n 
variable digits use 


x (ap) 31 Dividend 
yr ee), py(ap) Divisor 
a 


by e eecece Binary number that gives precision of y 
z(ap pz(ap) Quotient : 
pz(ap)  s_ ===--- Binary number that gives precision of z 
(sp) 31 Initially abs(x); successive remainders as 
algorithm progresses. 
stkz2(sp) d Temporarily holds the next d 
digits of quotient. ? 
stktl(sp) 31 sonporery because packed instructions 
don't allow overlapped operands 


OTS$DIV_PKSHORT 
1-801 


16-SEP-1 AX/VMS Macro Yet 


~3ee=18be P9EFEi33 HONSHTS Bact oO mars 
1 Holds abs(y 


OoOoOSClOSCOSOSCOOSCOSOSOSOSOSOSOSOOOCOOSOOOOOSOOOOOOOOSoOO 


a ee ee ca ee a ee ee ee ee a a ee a ed ed ad 
PUPP BS BS BS ES BE BE AWN nn nononononononunn ] 
NOUS UN OOD NOUS WN SDODNAU EWN "OO ONOAOU EWN OOONOM 


DOOCOCOSCSOOSOOOSOSOOSSOSOSOSOSCOSOOOOCOSOSOOOCOSOOOOOOOOOSSOOOOO 


DOOCSOSSOSS & FW SDooooo 
SOOOCOOOOOO COCO O OOOO OOO OSSSOSOSCO SOOO CCOOC COCO SCOOOOOOOOOCOOOOOOOO 


0 
0 
0 
0 
0 
Q 
0 
0 
0 
0 
0 
0 
Q 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


Soooooooo 

a 
PDAAOOOOUM 
QAuUilwr—Owooe 


So 


register usage: 


r10 z(a 


optimization notes: 


1) 
¢ Optimized for 


4 


a 


Soffset 0,,<- 
<,16 


~16>,- 
<stkz2,16>,- 
<stkt1,16>,- 
{6>,- 


<stksign, 1>,- 
<sthien,O>,< 


stkz2: 
stktl: 
stky: 
stksign: 
stklen: 


: parameter offsets 
; Soffset 4,,<- 


<x>,- 
<y>,> 


<py>,- 
<2o,- 
<p2>,- 


the sign of the 


ré6 a = additional digits of precision required beyond prec(x) 
r? stky(sp) = address of divisor 

r8 py(ap) = precision of divisor 

r9 r = number of additional digits of the quotient 
that are to be found for next step 


p) 
rij d = 31 = prec(y) = max. no of digits obtained each iteration 


Optimized for ey 5 not space. 


> Ud. 
for register to register operations are the same 
for byte operations and longword operations. 
Many packed instruction sequences were timed. Do not change 
unless actual tests are made to determine relative speed. 
Tests were made on 11/780 and Comet. 


Assumes speed 


stack offsets for work area 


uotient, 2 bits 
zlength of work area 


dividend by reference 


, ec(y) by valu 
zprec(z) by value 


OTSS$DIV_PKSHORT 16-SEP-1984 00:31: AX/VMS Macro v04-00 
1-001 aes 90:91:3) LIBRTL.SRCJOTSPKDIVS.MAR; 1 
167 <consta>,=- za by value 
168 <constd>,- 3d by value 
169 > 3 
008 fs 
00 ys 
$1 : 
014 pz: 
018 consta 
01C constd: 
000 170; 
8388 4 ; psect declarations 
00000000 158 : -psect _ots$code pic, usr, con, rel, lcl, shr, - 
0000 174 exe, rd, nowrt, long 
0000 175; 
0000 176; constant data area 
4 177 ; 
0c 444 178 zero: «packed +0 ;local packed decimal constant zero 
49 139 : local symbol definitions 
0000000F 48 136 bytes_to_sign=15 zbytes to sign for fixed decimal 31 
0001 184: 
0001 185 ; 
Boo) 186 ; 
CFFC 8053 is sentry ots$div_pkshort,*M<iv.dv.r2,r3.r4.r5,r6,r7,.r8,°9,r10,r11> 
0003 189 :initialize registers and temporaries 
SE BF AE 9E 0003 191° movab -stklen(sp),sp smake room for temporaries 
5A 10 AC DO 0007 136 mov | z(ap).r zsave address of quotient 
57 30 AE 9€ QOO00B 19 movab = stky(sp),r7 saddress of divisor 
58 OC AC DO QOOF 194 mov py(ap).r precision of divisor 
5B 1C AC 00 Q013 195 mov constd(ap),r11 id = 31 = prec(y 
56 18 AC DO Q017 196 movl consta(ap),ré6 za =_scale(z) + scale(y) - scale(x) 
1B 197 3 - 31 + prec(x) 
40 AE 9% aie 198 clrb gitslon (eg) iclear sign flag 
6—€ 04 8C 1F 34 1 199 movp #31,ax(ap), (sp) smove x, set gone. code 
26 14 OO 00 bgotr 50$ zbranch if x> 
1€ 19 ; 01 blss 40$ sbranch if x<0 
b Riise 
08 BC 58 DS AF O00 37 7 .” cmpp4 #0,zero,r8,ay(ap) rset condition code 
OA 13 s] beq sbranch if divide by 0 
14 aC 00 CBAF 00 00 8 ashp #0,#0,zero,#0.pz(ap),(r10) ;z = 0 
6 $08 ret 
BF AF 00 ce wT . 00 7 e 9 30S: divp #0,zero,#0,zero,pz(ap),(r10) ;cause divide by 0 
04 rt: 19 ret 3 
i I ix not 0, determine sign of x 
4 13; 
4 14 dos: 


OTS$D1V_PKSHORT 16*SEP=1984 00:31:31 VAX/VMS Macro V04-00 Page 5 
1801 gS EE 382 90:94:33 YTBRTL. SREJOTSPKDIVS MAR: 1 . (1), 
40 AE 88 incb stksign(sp) sset low order bit 
OF AE 9 decb bytes_to_sign(sp) ix <0 so make it positive 
idetermine sign of y 
zy may be 0 at this point 
: code optimized for y>0 
67 08 BC 5B 34 508: 8 ay(ap),(r7) into t 
mov r ap),(r move nto temporar 
oc 18 boeg «| OOk Sbranch if y> 20” 
40 A 3 incb stksign(sp) sset neg indicator 
58 AG AF 00 08 BC ? 2 subp6 = r8,ay ap),#0,zero,r8,(r7) sconvert to positive 


a, = & ww cmpp4 = #31, (sp), r8,(r7) px<y? 
2. @ blss 95$ sbranch if x<y, i.e. shift of d is o.k. 
iy <x 
6A 14 AC Fc; a > sf : divp 8. (r 7) #51. (sp) pz (ae) (r10) ;2z=x/y 
1F 6A 14 AC 5 38 5 mulp f8,(r7),pz(ap),(rl0) ,#31,stktlisp) :tl=(x/y) ey 
14 11 brb 110$ ; 


SSID DDD DDS IT 


ROD NEL HE EME BANU ODDO £9 MIWIIWIMMMMmMovmoowmwawawawawmwouw 


ax <y 
14 AC 00 83 AF O00 O00 F8 éss: ashp #0,#0,zero,#0,pz(ap),(r10) ;clear quotient 

; brb 115$ : 

;start of multiple precision divide 

100s: mulp r8,(r7),R11,stkz2(sp) ,#31,stktl(sp) ;tl=y*z2 


110$: pen #31,stktl(sp),#31,(sp) ;x=x-t1 
beq 150$ sbranch if remainder = 0 


7 Wee BD oF 3S 2 


:determine r, the number of the next low order digits to obtain 


DOOOOCOOOOOSOSOOSOSOSOOOOOOOOOOOOOOOOOOOOOOOOOO 


SOCCSCOCCOOOOCSoOSOOSoOoOooO 
OWOOOOOOO OCWDWWOMODdDcMcd ~~~ 


| 
| 
60$ 
start of divide proper; setup 


POPOPOPOPONONYD «= FOROROPOROPORORORORONONY «= FOROPORORY) «so ROROPOPORD «oo POPoronoronofroronorfnony:)o POnoTororonornornornorororyn» 
SOOOAAOAOAQ VNIWMMvmunuvmel es LLEEL LEEW WWW fporonronononon] ooo 
QuilWwn—O OOnNOuslWwN—OOW NOUFW "OOD NOULSWR—OOON CUS WNR)—OOO~aCWu 


59 5B 00 115$: movil ri1,r9 rr=d 
58 9 D1 cmpl r¢ rii 3a>d? 
0 14 bgtr 136$ ibranch if Larger 
59 28 4 movl 8.09 -r=a | 
20 AE 1F 00 GF iF 59 F 1308: ashp  £9,#31,(sp),#0,#31,stktl(sp) ;shift x left by r | 
6— 20 AE 1F 6434 A7 movp #31, stkel (sp) (sp) ;cOpy back into x 
14 AC 00 6A 16 ato 39 8 aC ashp 9, pz(ap),(r16),#0,pz(ap),stktl(sp) ;shift z left by r 
6A 20 AE 144A 4 6 mov (ap),stktl(sp),(r10) scopy back into z 
10 AE 5B CE 5 3 OBC divp Fa. (r7$ 831, (sp) rit astks (ap) szata)en/y | 
6A 14 AC 10 AE 38 C4 addp4 Ul estkz2(sp) .p2lap), (rt ) s2=2z42z | 
56 C cB subl2 r9,r jazacr | 
oS. F C bneq 106$ sbranch if more 
16 40 AE F8 oO bibs stksign(sp),155$ tbranch if quotient <0 | 
re ° 


TSS$D1V_PKSHORT 
-001 


0D 67 ; 
38 rt ; remainder = 0 
D 5 1508: 
14 AC 00 6A 16 AC 56 FB 0005 2371 ashp 
20 AE DD 
Ob 40 AE €8 OODF 2 blbs 
6A 20 AE 14 AC 4 E3 7 movp 
na 
EA A: ;quotient <0 
20 AE 6A AC 34 DEA 8 1558: movp 
FFO7 CF 00 20 AE 14 AC 23 dor 80 "  subp6 
6A 16 AC OF 
04 OOF C Hf ret 
OOFD 88 end 


of 
r6,pz(ap),(r10),#0,pz(ap), 


et pd get 160$ 
pz(ap),st ti(sp),(r10) 


pz(ap),(r10),stktl(sp) 


abs 9 00:31:31 woe Macro v04-00 
P-1984 11:15:22 CLIBRTL.SRCJOTSPKDIVS.MAR; 1 


emainder = 0 
stktl(sp) ;account for scale 


sbranch if quotient < 0 
;copy back into quotient 


scopy quotient into temp 


enter if tl holds quotient 


sen 
pz(ap),stktl(sp),#0,zero,pz(ap),(ri0) ;make z negative 


3 source lines were read in Pass 1, producin - ect records in Pass 2. 
pages of virtual memory were used to define AE 


BSE Sowers nr ee swe ne sma eeres } 


OTS$DIV_PKSHORT 16-SEP- AX/VMS Macro V04-00 Page 7 | 
Symbol fable oreEpm bbe ii ssas LO NEN TS eae s eet Cs mans1 729 ay 
BYTES TO_SIGN = Q000000F 
CONSTA ~ 09000018 
CONSTD 00 Ooi 
DIR... E i 5 1 
OTSSDIV_PKSHORT 00000001 RG 02 
PY 000000C 
464 
STKLEN 0000041 | 
STKSIGN 0000040 
STKT1 BD b OSS 
STKY 000000 
STKZ2 aid 
Xx 0000004 
Y 00000008 | 
Zz 00000070 
ZERO 00000000 R 02 
foc emeoceanenenoeoee + 
! ;_Psect synopsis : 
PSECT name Allocation PSECI No. Attributes 
. ° 00000000 < 0.) 00 ¢ O.) NOPIC USR CON ABS LCL te NOEXE NORD NOWRT NOVEC BYTE 
SABSS 00000047 65.) 01 ¢ 1.) NOPIC USR CON ABS - LCL NOSHR EXE RD WRT NOVEC BYTE 
_OTSSCODE OOOOOOFD ¢ 253.) 02 ¢ 2.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
OOO EEE B OM ew remem ee + 
! ; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 39 00:80:00 -09 00:00:02.69 
Command processing 117 00:00:00.31 00:00:0 £6 
Pass 116 00:00:00.80 0:00:05.2 
Symbol table sort 83:83: 0.01 3:88: 0.01 
Pass 2 6 0:00:00.44 0:00: 5-28 
Symbol table output 00:00:00.02 0:00:00.02 
Psect synopsis output 00:00: fat 0:00:00.01 
Cross-reference output ! 00:00:00.00 0:00: 9-90 
Assembler run totals 34 00:00:01.64 0:00:15.18 
ihe working set Limit was 900 pages. 
123 byies (14 pages) of virtual memory were used to buffer the intermediate c sgd e. | 
sd were 10 pages of symbol table space oh ecas d to hold 17 non-local and 12 local symbols. 
| 
| 
| 


Macro Library name Macros defined 


~-$255$DUA28: CSYSLIBISTARLET.MLB; 2 


M 2 
VAXot1 Racto Run Statistics NOEEH TERS i FEiES EMGRTE. SRESBreRROIGs makes P29 ef, 
45 GETS were required to define 2 macros. 


There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPKDIVS/OBJ=0BJ$:OTSPKDIVS MSRC$:OTSPKDIVS/UPDATE=(ENH$:OTSPKDIVS) 


071 AH-BT13A-SE EQUIPMENT CORPORATION 
VAXZVMS V4.0 NTIAL AND PROPRIETARY 


